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Abstract. We exemplify the way the rigged Hilbert space deals with the Lippmann- 
' Schwinger equation by way of the spherical shell potential. We explicitly construct the 

04 ■ Lippmann-Schwinger bras and kets along with their energy representation, their time 

evolution and the rigged Hilbert spaces to which they belong. It will be concluded 
^ ' that the natural setting for the solutions of the Lippmann-Schwinger equation — and 

\Q , therefore for scattering theory — is the rigged Hilbert space rather than just the Hilbert 

space. 

cn 
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^ ■ PACS numbers: 03.65.-w, 02.30.Hq 

^' 

Qh, 1. Introduction 

+-> 

■ Scattering experiments play a central role in our understanding of the microscopic world. 
For example, Compton scattering of x rays by electrons is held as experimental evidence 
for the particle nature of the photon; much of what has been learned about the structure 

^ ■ of the nucleus, indeed even its discovery, was the results of scattering experiments; the 

■ analysis of scattering has yielded most of our present knowledge of elementary particle 
physics [1]. Scattering theory was developed in order to understand such experiments. 

The Lippmann-Schwinger equation is one of the cornerstones of scattering theory. 
It was introduced by Lippmann and Schwinger in 1950 [2], and it has become a 
commonplace in scattering theory [3-6]. It is written as 

where \E^) represent the "in" and "out" Lippmann-Schwinger kets, \E) represents an 
eigenket of the free Hamiltonian Hq, 

Ho\E) = E\E), (L2) 

and V is the potential. The Lippmann-Schwinger kets are, in particular, eigenvectors 
oiH: 

H\E^) = E\E^). (L3) 
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To the kets \E^), there correspond the bras (^-El, which satisfy 

{^E\ = {E\ + {^E\V I (1.4) 

The bras (^-E] are left eigenvectors of if, 

{^E\H = E{'^E\, (1.5) 
and the bras {E\ are left eigenvectors of Hq, 

{E\H = E{E\. (1.6) 

In spite of its wide acceptance, the Lippmann-Schwinger equation still lacks a proper 
mathematical setting. There are mainly two problems. First, the bras {^E\ and the kets 
l^^^*") are not square integrable (i.e., they lie outside the Hilbert space), and therefore 
they must be treated as distributions. And second, in resonance theory the bras {^E\ 
and the kets \E^) need to be analytically continued into the complex plane [7-14]. Both 
problems can be handled by means of the rigged Hilbert space (RHS). In this paper 
we shall focus on the first problem, namely on how to accommodate the Lippmann- 
Schwinger bras and kets for real, positive energies within the RHS setting. The second 
problem will be addressed in a sequel to this paper. 

The present paper is the continuation of Refs. [15-18], where we showed how the 
RHS deals with continuous spectrum. As we shall see, the methods of Refs. [15-18] can 
be straightforwardly applied to the Lippmann-Schwinger equation. 

We recall that a rigged Hilbert space is a triad of spaces (see Ref. [19] for a 
pedagogical introduction to the RHS) 

$C7^C$^ (1.7) 

such that 7"^ is a Hilbert space, $ is a dense subspace of H, and is the space of 
antilinear functionals over The space is called the antidual space of The 
space $ is the largest subspace of the Hilbert space such that first, it remains invariant 
under the action of the Hamiltonian, and second, the action of the Lippmann-Schwinger 
bras and kets is well defined on its elements. Associated with the RHS ()1.7j) . there is 
always another RHS, 

^CHC^', (1.8) 

where is called the dual space of $ and contains the linear functionals over As we 
shall show, the Lippmann-Schwinger bras {^E\ belong to whereas the Lippmann- 
Schwinger kets \E^) belong to 

It should be emphasized that our RHS does not entail an extension of the physical 
principles of scattering theory. Our RHS simply enables us to extract and process 
information contained in the Lippmann-Schwinger bras and kets. 

There have been previous attempts at describing the Lippmann-Schwinger equation 
within the RHS setting in Refs. [8-14] . Our results will be different from those of [8-14] 
in the following sense. First, Refs. [8-14] use formal manipulations based on the Hardy- 
space assumption, whereas our approach is based on an explicitly solvable example; 
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what is more, we do not make any assumption concerning Hardy functions, because 
Hardy functions seem unrelated to the solutions of the Lippmann-Schwinger equation. 
Second, Refs. [8-14] use two different RHSs (of Hardy functions), whereas we shall use 
only one RHS. Third, the Hardy-function assumption of Refs. [8-14] implies that the 
time evolution of the Lippmann-Schwinger bras and kets associated with real energies 
is given by a semigroup, whereas in this paper the time evolution of the Lippmann- 
Schwinger bras and kets is given by a group. And fourth, using Hardy functions as in 
Refs. [8-14] entails an extension of the physical principles of scattering theory, whereas 
in this paper we shall remain within standard quantum scattering theory. 

The paper is organized as follows: In Sec. El we recall some of the basics of 
scattering theory. In Sec. 01 we obtain the Lippmann-Schwinger eigenfunctions in the 
radial, position representation for the spherical shell potential and for zero angular 
momentum. In Sec. (H we recall the expressions for the domain, the spectrum and the 
Green function of the Hamiltonian. In Sec. we obtain the unitary operators U±, 
the eigenfunction expansions and the direct integral decompositions generated by the 
Lippmann-Schwinger eigenfunctions. In Sec. (HI we construct the Lippmann-Schwinger 
bras and kets from the corresponding eigenfunctions, and we obtain their most relevant 
properties. In Sec. [71 we calculate the time evolution of the Lippmann-Schwinger bras 
and kets. In Sec. |H1 we justify some standard results of scattering theory within the 
RHS setting. In Sec. O we state our conclusions. 

2. A reminder of scattering theory 

A scattering process can be pictorially summarized as in Fig. ^ Loosely speaking, we 
send a beam of prepared initial "in" states toward the potential. After the collision 
takes place, becomes ^9°"*. We then measure the probability to find a final "out" 
state ip°^^ . The amplitude of this probability is given by 



where 5* is the S matrix. The state <y9™ is determined by initial conditions through a 
preparation procedure (e.g., an accelerator). The state ^9°"* = Sip^'^ is determined by 
those initial conditions through and by the scattering process through the S matrix. 
The state -j/^""* is determined by final conditions through a registration procedure (e.g., 
a detector). 

The initial "in" state ip^^ and the final "out" state ip"^^ are asymptotic forms of the 
so-called "in" state y?"*" and "out" state ip~ in the remote past and in the distant future, 
respectively. In terms of these, the probability amplitude (j2.Hl reads 



The asymptotic states (p^^ and ip"^^ are related to the "exact" states y^"*" and ip by the 
M0ller operators. 



(2.1) 



(2.2) 



in 




(2.3a) 



out 



(2M) 
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It is important to keep in mind that the states <y9™, v?""* and ip"^^ (or if^ and are 
not the central object of quantum scattering. The central object of quantum scattering 
is the probability amplitude ()2.Hl - ()2.2j) . which measures the overlap between the "in" 
(i.e., prepared) and the "out" (i.e., detected) states. This is in distinction to classical 
scattering, where the states ^9™, and 99°"* do actually represent a physical wave, and 
where we actually are able to observe 93°"* rather than the overlap of <^°"* with some 

^out 

It is customary to split up the (total) Hamiltonian H into the free Hamiltonian Hq 
and the potential V , 

H = Ho + V. (2.4) 

The potential V is interpreted as the interaction between the components of the initial 
prepared states, for instance, the interaction between the in-going beam and the target. 
The states ^9™ and evolve under the influence of the free Hamiltonian Hq, whereas 
the states and ip~ evolve under the influence of the (total) Hamiltonian H. 

The dynamics of a scattering system is therefore governed by the Schrodinger 
equation subject to boundary conditions that specify what is "in" and what is "out." 
The Lippmann-Schwinger equation for the "in" \E~^) and "out" l-E"") kets has those 
"in" and "out" boundary conditions built into the ±ie, since Eq. fll.l|l is equivalent to 
the time-independent Schrodinger equation ()1.3|) subject to those "in" {+ie) and "out" 
(— ie) boundary conditions. In the position representation, the +ie prescription yields 
the following asymptotic behavior: 

(x|E+) ,e''^ + f{e,<p) — , (2.5) 

r^oo r 

where x = (x, y, z) = (r, 6*, 0) are the position coordinates, k is the wave number and 
f{9, (p) is the so-called scattering amplitude. Thus, far away from the potential region, 
(x|i?+) is a linear combination of a plane wave (which originates from the free part in 
Eq. (jl.ll) ) and an outgoing spherical wave multiplied by the scattering amplitude (which 
originate from the second term on the right-hand side of Eq. (jl.ip . and which account 
for the effect of the potential on the incoming beam). The —ie prescription leads to the 
following asymptotic behavior: 

(x|E-) > e'^^ + 7(M) ^ ; (2.6) 

r^oo T 

thus, far away from the potential region, (x|ii^~) is a combination of a plane wave and 
an incoming spherical wave multiplied by the complex conjugate of f{9, (p). 

Formally, the "in" state (p~^ and the "out" state ip~ can be expanded in terms of 
the Lippmann-Schwinger bras and kets as follows: 

POO 

cp+= / dE\E^){+E\v+), (2.7) 
Jo 

POO 

ip~ = dE\E~){~E\i)-) . (2.8) 
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The "in" ("out") states can be also expanded in terms of the "out" ("in") Lippmann- 
Schwinger bras and kets: 



/ dE\E-){~E\ip+), (2.9) 

^0 

POO 

^- = I dE\E+){'^E\i>-) . (2.10) 

By combining Eqs. ()2.7p - ()2.8|) with the relation 

{-E'\E+) = S{E)5{E - E') , (2.11) 
we can express the matrix element ()2.2j) as 

POO 

(^-,¥,+)=/ dE{r\E-)S{E){-^E\^^). (2.12) 

Formally as well, the initial and final states can be expanded by the bras and kets 
of the free Hamiltonian: 

POO 

/ dE\E){E\^^), (2.13) 

POO 

= / dE|E)(E|V^°"*), (2.14) 
and the probability amplitude ()2.H) can be written as 

POO 

(^°"\V)=/ dE(V^°"*|E)S(E)(E|(/?'°), (2.15) 

where we have used Eqs. ()2.13j) - ()2.14j) and the relation 

{E'\S\E) = S{E)5{E-E'). (2.16) 

As we shall see in Sec. the formal expressions ()2.7j) - ()2.1(i|l acquire meaning within the 
RHS. 

3. Solving the Lippmann-Schwinger equation 

For the spherical shell potential, the Lippmann-Schwinger equation can be solved 
explicitly. For the sake of focus, we shall restrict ourselves to angular momentum / = 0. 

3.1. The radial Lippmann-Schwinger equation 

The expression for the spherical shell potential is given by 

{0 < r < a 
Vo a<r<b (3.1) 
6 < r < oo , 

where Vq is a positive number that determines the strength of the potential, and a 
and b determine the positions in between which the potential is non-zero. Because 
the potential ()3.H1 is spherically symmetric, we shall work in the radial position 
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representation. In this representation and for / = 0, the free Hamiltonian Hq acts 
as the formal differential operator Hq, 

Hofir) = h,fir) = -^§-Jir) , (3.2) 
V acts as multiplication by the rectangular barrier potential V{r), 

{Vf){r) = V{r)f{r), (3.3) 
and the total Hamiltonian H acts as the formal differential operator h, 

Hf{r) = hf{r) = (-|^^ + Vir)^ f{r) . (3.4) 

In the radial representation, the Lippmann-Schwinger equation becomes 

(r|E±) = {r\E) + {r\—^-—-V\E^) . (3.5) 

In Eq. ()3.5|) . the {r\E) are eigenfunctions of the formal differential operator Hq, 

ho{r\E) = E{r\E) , (3.6) 

whereas the {r\E^) are eigenfunctions of the formal differential operator h, 

h{r\E^) = E{r\E^) , (3.7) 

subject to the boundary conditions specified by Eqs. ()3.9aj) - fl3.9e|) . 

3.2. The solutions to the radial Lippmann-Schwinger equation 

The procedure to solve Eq. ()3.5|) is well known [5,6]. Since Eq. fl3.5p is an integral 
equation, it is equivalent to a differential equation subject to the boundary conditions 
that are built into that integral equation. In our case, for / = 0, Eq. ()3.5|) is equivalent 
to the Schrodinger differential equation, 

+ Vir)] {r\E^) = E{r\E^) , (3.8) 



2m dr^ 

subject to the following boundary conditions: 

{r\E^) = 0, (3.9a) 
{r\E^) is continuous at r = a,b, (3.9b) 

-^{r\E^) is continuous at r = a,6, (3.9c) 
dr 

{r\E+) ~ e-'''' - S{E) e'^"^ as r ^ oo , {3.9d) 



{r\E-) ~ e"^'' - S{E) e'^^'' as r ^ oo , (3.9e) 



where 



k = f-^E (3.10) 

is the wave number, and S{E) is the S matrix in the energy representation. We recall 
that the boundary conditions ()3.9(jj) and (j3.9ep originate from the -\-ie and from the —ie 
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conditions of Eq. (j3.5p . respectively. The boundary condition ()3.9(jj) means that far away 
from the potential region, {r\E^) is a combination of an incoming spherical wave and 
an outgoing spherical wave multiplied by the S matrix. The boundary condition ()3.9ej) 
means that far away from the potential region, {r\E~) is a combination of an outgoing 
spherical wave and an incoming spherical wave multiplied by the complex conjugate 
of the 5* matrix. The asymptotic behaviors ()3.9d|) and ()3.9ej) are the / = 0, radial 
counterparts of the asymptotic behaviors ()2.5|) and ()2.6|) . 

If we insert expression ()3.1|1 for the potential into Eq. ()3.8|1 . and solve Eq. ()3.8j) 
subject to the boundary conditions ()3.9a)) - ()3.9ep . we obtain the "in" and "out" 
eigenf unctions, 

{r\E^) ^ xHr; E) = iV(E) , Ee[0,oo), (3.11) 

where N{E) is a delta-normalization factor. 




x{r; E) is the so-called regular solution of Eq. (|3.8p 

x{r]E) = 



sin(y r) < r < a 

J,{E)d^^^'- + ME)e-'^^^' a<r<b (3-13) 
J3(E)eV^^" + J4(E)e"'v^" 6<r<oo, 
and J±{E) are the Jost functions, 

J+{E) = -2iJi{E) , (3.14a) 

J_{E)=2iJ^{E). (3.146) 

The explicit expressions for J\-J^ can be obtained by matching the values of xij] E) 
and of its derivative at the discontinuities of the potential (see Eqs. flB.l|) - f)B.4|) in 
[Appendix B| ). In terms of the Jost functions, the S matrix is given by [5,6] 

S(E) = |-|||, i?G[0,oo). (3.15) 

From Eq. ()3.11|) it follows that the "in" and "out" Lippmann-Schwinger eigenfunctions 
are proportional to each other, 

X^ir;E)=S{E)x-{r;E). (3.16) 
3.3. The "left" Lippmann-Schwinger eigenfunctions 

When we write Eq. ()1.4|1 in the radial position representation, we obtain an integral 
equation for the "left" Lippmann-Schwinger eigenfunctions, 

{^E\r) = {E\r) + {^E\V !r ^ J r) . (3.17) 



+ V{r)] {^E\r) = E{^E\r) , (3.18) 
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Solving Eq. (|3.17p is analogous to solving Eq. (j3.5|) . For / = 0, Eq. (|3.17|) is equivalent 
to the Schrodinger differential equation, 

2m dr^ 

subject to the following boundary conditions: 

{'^E\0) = 0, (3.19a) 

(^i?|r) is continuous at r = a,6, (3.19&) 

— C^E\r) is continuous at r = a,6, (3.19c) 
dr 

{+E\r) ~ e''^" - SCE) e-''^" as r ^ oo , (3.19(i) 

{-E\r) r^e-"""" - S{E)e"'' as r oo . (3.19e) 

Note that the boundary conditions ()3.19dj) and (j3.19ep originate from the —ie and from 
the +ie conditions of Eq. ()3.17p . respectively. 

By inserting Eq. 1)3.11) into Eq. ()3.18p . and by solving Eq. ()3.18|) subject to the 
boundary conditions ()3.19a|) - fl3.19e|) . we obtain that the "left" Lippmann-Schwinger 
eigenfunctions {^E\r) are the complex conjugates of the "right" Lippmann-Schwinger 
eigenfunctions {r\E^), 



{^E\r) = x^{r;E). (3.20) 



4. Domain, spectrum and resolvent 

In this section, we include the expressions for the domain, the spectrum and the resolvent 
of H. Such expressions were obtained in Ref. [14], and are included in this section 
because we shall need them in subsequent sections. 

The formal differential operator ()3.4|) has an infinite number of self-adjoint 
extensions [20]. These self-adjoint extensions are characterized by the following 
boundary conditions [20]: 

/(O) + «/'(0) = 0, -oo<«<oo. (4.1) 
Among all these, the boundary condition needed in scattering theory is 

/(0) = 0. (4.2) 

The boundary condition ()4.2j) selects the following domain for the Hamiltonian: 

V{H) = {/(r) I /(r), hf{r) G ^^([0, oo), dr), /(r) G AC'[0, oo), /(O) = O} , (4.3) 

where ylC^[0,oo) denotes the space of functions whose first derivative is absolutely 
continuous. This domain induces a self-adjoint operator H, 

W)(r)=(-|^^ + nr))/(r), feV{H). (4.4) 

The spectrum of H, Sp{H), was shown in Ref. [14] to be [0, oo). In Ref. [14], we 
also calculated the Green function G{r, s; E) of H for different regions of the complex 
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energy plane. In all our calculations, we used the following branch for the square root 
function: 

: {E eC\ - TT < aTg{E) < tt} i — > {E eC\ -n/2< aTg{E) < n/2} . (4.5) 
The branch ()4.5p grants the following relation: 



E =VE, EeC, (4.6) 

which in general does not hold for other branches of the square root function. 

Either by applying Theorem 1 of [Appendix A[ or by borrowing the results from 
Ref. [14], one can see that in the first quadrant of the energy plane, the Green function 
can be written as 

r -nNiE)x+{r;E)f+{s-E) r<s 
G{r,s;E)={ Re(E) > , Im(E) > , (4.7) 

[ -TxN{E)x-^{s-E)f+{r-E) r>s 

whereas in the fourth quadrant its expression reads 

r -nN{E)x-ir;E)f-{s;E) r<s 
G{r,s;E)={ Re(E) > , Im(E) < . (4.8) 

[ -TxN{E)x-{s-E)f-{r-E) r>s 

In Eqs. (gUjl-dlSI), X^{r\ E) are given by Eq. (jSHH), N{E) is given by Eq. (jTT^ . and 
the eigenfunctions /^(r; E) are given by 



^±(^)eV^^^' + ^±(E)e~V7^^'^ 0<r<a 
fHr; E) = { Af{E)l^^^" + A^(E)e-V¥(^^ a<r<b (4.9) 
e V ft-i < r < oo . 

The coefficients Af-A^ can be obtained by matching the values of /^(r; E) and of their 
derivatives at the discontinuities of the potential (see Eqs. ()B.6|1 - (IB.9|1 in [Appendix B ). 



5. The operators U± and Uq 

In this section, we obtain the Fourier-like transforms U±, the eigenfunction expansions 
and the direct integral decompositions generated by the "in" and "out" eigenfunctions. 
The operators U±. will be obtained by applying the Sturm-Liouville theory (see Ref. [20] 
and Appendix A ) . We shall end this section by recalling the expression for the Fourier- 



like transform Uq associated with the free Hamiltonian. 

5.1. The "in" unitary operator U+ 



The procedure to obtain f/+ consists of applying Theorems 2, 3 and 4 of Appendix A 
In order to be able to apply Theorem 4, we choose the following basis for the space of 
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solutions of ha = Ea that is continuous on (0, oo) x (0, oo) and analytically dependent 
on in a neighborhood of (0, oo): 

ai(r;E) = x+(r;^), (5.1a) 



cos[ ^/^Er) 



(T2{r-E) = { 



+ C2{E)e 



-h TT(.E-Vo) r 



■ /2rr 



^'^ + C4(i?)e-V^^^ 



< r < a 
a<r<h (5-1^) 
6 < r < oo . 



The functions C1-C4 can be obtained by matching the values of (J2{r] E) and its first 
derivative at the discontinuities of the potential (see Eqs. f)B.10|) - ()B.13|) in [Appendix B ). 
Equations (gSl) and (IKTall - fl5.16jl lead to 



2iME)C^{E) , Ja{E) 

f+{r; E) = —————— ai(r; E) + — — -(T2(r; E) 



and to 



where 



rir;E) 



N{E) W{E) 

2iME)C,{E) 
N{E)W{E) 



W{E) 

a,{r;E)-^j^a2{r;E) 



W{E) 



W{E) = J^{E)C^{E) - ME)C^{E) . 
After substituting Eq. ()5.2j) into Eq. ()4.7j) we arrive at 



G{r,s-E) = -71 



N{E) 
W{E) 



2ij4(E)C4iE) 



>-ai{r; E) + ME)a2{r- E)J ai(s; E) , 

Re(^) > , \m{E) > , r>s. 
After substituting Eq. into Eq. ()4.8j) we arrive at 



G{r,s]E) = -7T 



NjE) Ji{E) 
W{E) ME) 



E) + UEy2{r- E)] a,{s- E) , 



Re(E) > , \m{E) < , r> s 



where we have used the relation 

ME) 



X {r;E) 
Because by Eq. fl4.(i|l 



x+{s;E) 

Eq. (jniSD leads to 



ME) 
ME) 



X^{r;E). 



X-'is;E) 



G{r, s- E) = 27ri^^§^ a,{r- E)a,{s; E) + ^^(r; E)a^is; E) 



Re{E) > , lm{E) > , r>s, 



and Eq. (j5.6j) leads to 



G(r, E) = 27ri^^^§p ai(r; E)a,{s- E) + TrMfP ^^(r; E)cr,{s- E) , 

Re{E) > , lm{E) < , r > s . 



(5.2) 

(5.3) 
(5.4) 

(5.5) 



(5.6) 



(5.7) 



(5. 



(5.9) 



(5.10) 
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On the other hand, the Green function can be written in terms of the basis {ai, (T2} of 
Eqs. flKTaj) - (l5T^ as (see Eq. (TO!) in Theorem 4) 

2 

G{r, s; E) = ^ e+{E)a,{r- E)a,{s; E) , r>s. (5.11) 

By comparing ()5.11|1 to we obtain 

Ot^E) = NiE')}^(^E) . , ME) > , Im(i?) > . (5.12) 

\ ^ W(E) U / 

By comparing (j5.1ip to (j5.1(Jj) . we obtain 

/ 2T,i ME)C,{E) Q \ 

= ^(4^3^],) , Re{E) > , lm{E) < . (5.13) 

\ ^ W{E) U / 

From Eqs. (j5.12p and (j5.13j) it follows that the measures pi2, P21 and P22 of Eq. ()A.10|1 

are zero, and that pn is given by 

n ((F F))- ^ r 2^- ME)Cs{E) - ME)C,{E) 

-E2 

dE = E2-E^, (5.14) 

'El 

where we have used (|5.4|) in the second step. Therefore, the measure pn is just the 
Lebesgue measure, which means, in particular, that the eigenfunctions x'^i'^^'^E) are 
5-normalized (see also Sec. 2.9 of Ref. [16]). 

By Theorem 2 of [Appendix A[ there is a unitary operator U+ that transforms from 
the position representation into the energy representation, 

[/+ : L2([0,cx)),dr) 1 — > L\[0,oo),dE) 

(5.15) 

/(r) ^ /+(E) = (f/+/)(E) = /;°dr/(r)x+(r;E), 

where f+{E) denotes the energy representation of the function /(r) when obtained by 
way of f/+. The action of f/+ on the domain T>{H) is given by 

V{H) = U+V{H) = !^UiE) e L\[0, 00), dE)\ ^ dE |E/+(E)|2 < ooj . (5.16) 

It can be easily checked that the operator f/+ diagonalizes H in the sense that 
H = U+HUl^^ acts as the operator multiplication by E, 

H -.ViH) C L'^{[0,oo),dE) 1 — > {[0 , 00) , dE) 

U ^ HU{E)=EU{E). ^^'^^^ 

The inverse of U+ is given by Eq. dH)): 

poo 

f(r) = U-'Uir)= dEU{E)x+ir;E), ^E) E L\[0,oo),dE) . (5.18) 
^0 

The operator U^'^ transforms from the energy representation back into the position 
representation. Expressions ()5.15|) and ()5.18|) provide the eigenfunction expansions of 
any square integrable function in terms of the "in" eigensolutions. 
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The operator U^, and therefore the "in" eigenfunctions x'^i'^^'iE), entails a direct 
integral decomposition of the Hilbert space in a straightforward manner: 

n I — > u+n = n = ® j^^^^^ n{E)dE 
f u^f^\u{E)}, feV{H)J4E)eH{E). 



(5.19) 



In this equation, the Hilbert spaces H, H and Hi^E) are respectively realized by 
L2([0, oo), dr), L2([0, oo), dE) and C. 



5.2. The "out" unitary operator f/_ 

The construction of f/_ follows the same procedure as the construction of f/+. The 
functions 



ai{r;E) = x {r; E) 



(5.20a) 



cos 



^E r) 



< r < a 



(T2{r;E) = < 



Ci(E)eV^-^o) ^ + C,(E) ' a<r<b 
C3(^)eV^^'- + C4(E)e"V^^'' fe<r<oo 



(5.206) 



form another basis for the space of solutions of ha = Ea that is continuous on 
(0, oo) X (0, oo) and analytically dependent on in a neighborhood of (0, oo). Therefore, 
we are allowed to apply Theorem 4 of [Appendix A Equations ()4.9|1 and ()5.20aj) - (j5.20?)j) 
lead to 



/*(r;B) 



N{E)W{E) 



W{E) 



and to 



2iME)Cs{E) JsjE) 
f {r;E)= ai{r; E) - a2{r; E) . 



N{E)W{E) ' ' ' W{E) 
By substituting Eqs. ()5.21|) and ()5.7j) into Eq. ()4.7j) . we obtain 



G{r,s]E) = nN{E) 



ME) 



Re{E) > , Im(^) > , r>s. 
By substituting Eq. (j5.22p into Eq. (j4.8p . we obtain 



G{r,s;E) = -nN{E) 



Because by Eq. fl4.6|) 



2ij3(g)C3(g) 

NiE)W(E) 



a,ir; E) - §^a2ir; E)\ a,is; E) 
Re{E) > , lm{E) < , r > s 



X-{s;E) = -^^X'{s;E) 



ME)' 



Eq. (jK:^ leads to 



G(r, E) = 27ri^2^fJ|pai(r; E)ai(.; E) - vr^^^f|pa2(r; E)a^i.s; E) , 

Re(E) > , lm{E) > , r > s 



(5.21) 



(5.22) 



(5.23) 



(5.24) 



(5.25) 



(5.26) 
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and Eq. (IH^ leads to 



G(r, s; E) = 27ri^^^^§^a,{r; E)a,{s; E) - n^^^§^a,{r- E)a^{s; E) , 

Re(E) > , \m{E) < , r> s. 

On the other hand, by way of Eq. ()A.9|) . we can write the Green function in terms of 
the basis {ai, 02} of Eqs. (l5.2()ajl - flO(I^ as 

2 

G(r, s- E) = Otj{E)a,{r- E)a,{s- E) , r>s. (5.28) 

By comparing (j5.28p to (j5.26|) . we obtain 

/ ^TTiMMdEl 

G^E) = NjElSiE) . 1 , ME) > , lm{E) > . (5.29) 

V w{E) 

By comparing ()5.28|) to ()5.27p . we obtain 

/ 27riiM^1^3M 

OUE)=[ _ rmiiiEX . h Re{E) > , lm{E) < . (5.30) 

\ ^ WiE) U 

From Eqs. ()5.29p . ()5.30|) and ()A.10|) . it follows that the measures pi2, P21 and P22 in 
Theorem 4 of [Appendix A are zero, and that pu is given by the Lebesgue measure, 

pn((Ei,E2))= / dE = E2-E,. (5.31) 

JEi 

This means, in particular, that the x~{^'i E) are 5-normalized. 

By Theorem 2 of [Appendix A[ there is a unitary operator f/_ that transforms from 
the position into the energy representation, 

/■oo 

f_{E) = {U^f){E)= drf{r)x-{r;E), f{r)eL\[0,oo),dr), (5.32) 

Jo 

where f-{E) denotes the energy representation of /(r) when obtained by way of f/_. 
The inverse of U- is given by Eq. (jSH): 

/>oo 

f{r) = {Uz't){r)= dEf4E)x-{r;E), /l(i?) G ^^([0, 00), di?) . (5.33) 

Jo 

Likewise the operator UZ^ transforms from the energy representation into 

the position representation. Likewise U+, U- carries the domain 'D{H) onto the 
domain ()5.16[) . and H = U^HUZ^ acts as the operator multiplication by E. As well, 
expressions fl5.32|l and ()5.33|1 provide the expansions of any square integrable function in 
terms of the "out" eigenf unctions, and a direct integral decomposition similar to ()5.19|) . 

5.3. The "free" unitary operator Uq 

The operator Uq was constructed in Refs. [14,17]. Because we shall need Uq in order to 
construct the M0ller operators, in this subsection we recall the expression for Uq. 
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The regular (i.e., vanishing at r = 0), (5-normahzed eigensolution of the differential 
operator ()3.2|) is given by 



Xo{r]E) = N{E)sin\ \j'^E r\ , < r < oo . (5.34) 

This eigensolution can be used to construct the unitary operator 

t/o : L2([0,oo),dr) i — > L2([0, oo), dE) 

^ (5.35) 
/ Uof^fo 

that transforms from the position into the energy representation. The action of Uq can 
be written as an integral operator: 



foiE) = UofiE) = / dr/(r)xo(r; E) , /(r) G L^([0, oo), dr) , (5.36) 







where fo{E) denotes the energy representation of /(r) when obtained by way of Uq. 
The inverse of Uq can also be written as an integral operator: 

POO 

fir) = U,'foir)= dEfoiE)xoir;E), /o(E) G ^^([0, oo), dE) . (5.37) 
Jo 

Expressions fl5.36p and ()5.37|) provide the expansions of any square integrable function 
in terms of the "free" eigenfunctions. 

Note that, since the Lippmann-Schwinger eigenfunctions tend to the "free" 
eigenfunctions when the potential vanishes, 

\im x^ir;E) = Xoir;E) , (5.38a) 

Vo— >0 

the operators U±. tend to Uq when the potential vanishes, 

limt/± = f/o. (5.386) 

Vb^O 

6. Construction of the "in" and "out" bras and kets 

The solutions to the Lippmann-Schwinger equation are eigenvectors of the Hamiltonian 
whose eigenvalues lie in the continuous spectrum. As explained throughly in Refs. [14, 
19], eigenvectors whose eigenvalues lie in the continuous spectrum must be treated as 
distributions by way of the rigged Hilbert space. In this section, we construct the "in" 
and "out" bras and kets together with the rigged Hilbert spaces that accommodate 
them. As it turns out, the RHS constructed in Refs. [14-16] suffices for such purpose. 
The results of this section will be summarized by Proposition 1 at the end of this section. 

6.1. The "in" kets 

The definition of a ket is borrowed from the theory of distributions [21]. Given a function 
f{x) and a space of test functions the antilinear functional F associated with the 
function f{x) is an integral operator whose kernel is precisely f{x): 



F{ip)= / dxip{x)f{x). (6.1) 



The RHS approach to the Lippmann-Schwinger equation 



15 



The "bad behavior" of the distribution f{x) must be compensated by the "nice behavior" 
of the test function V5(x), so the integral ()6.ip makes sense. 

By using definition (jfi.lj) . we associate an "in" ket \E^) with the "in" eigenfunction 
X+(r; E) for each E G [0, oo): 

\E+) : $ I — ^ C 



In Dirac's notation, the action of \E^) is written as 



(6.2a) 



{^^\E+) = dr ((^+|r)(r|E+) . {6.2b) 







Note that even though E) = {r\E~^) is also meaningful for complex energies, the 

energy in Eqs. (j6.2a|) - ()6.26|) runs only over Sp(if) = [0,cxd), because in this paper we 
restrict ourselves to bras and kets associated with energies that belong to the spectrum 
of the Hamiltonian. 

We now need to find the subspace $ on which definition ()6.2a|) makes sense. Besides 
making ()6.2aj) well defined, the space $ must also be invariant under the action of the 
observables of the system. The invariance of $ is a crucial property, since it entails 
finite expectation values and well-defined commutation relations, and since it allows us 
to apply observables on the elements of $ as many times as wished [19]. Since in this 
paper the only observable we are concerned with is the Hamiltonian, we shall simply 
require invariance under H. Thus, the space $ must satisfy the following conditions: 

• The space $ is invariant under the action of H. (6.3a) 

• The elements of $ are such that the integral in Eq. ()6.2aj) makes sense. {6.3b) 

In order to meet requirement ()6.3a| . the wave functions ^p^{r) must at least be in 
the maximal invariant subspace of H: 

oo 

1) = Pi V{H'') . (6.4) 

n=0 

As shown in Refs. [14, 16], the space V is given by 

V = {ipe L2([0, oo), dr) I /iXr) G ^^([0, oo), dr), (/?(")(a) = (/?(")(6) = 0, 

/iX0)=0,n = 0,l,2,...; (^(r) G C-([0, oo))} , ^ ' ^ 

where y?'-"^ denotes the nth. derivative of ^p. In order to meet requirement (j6.3&|l . the wave 
functions ^p~^{r) must behave well enough so the integral in Eq. ()6.2aj) is well defined and 
yields a continuous, antilinear functional. From the expression for x^('"j E), Eq. ()3.1H1 . 
one can see that the <f~^{r) have essentially to control purely imaginary exponentials. 
Therefore, the space $ constructed in Refs. [14,16] meets the requirements ()6.3aj) - ()6.35|) : 

G L2([0,oo),dr) I G y^\\n,m < oo , n, m = 0, 1, 2, . . .} , (6.6) 

where the || \\nm are given by 



y^'lkrn ■■= xl I dr \{l + r)-{l + H)"^ip+{r)f , n,m = 0,1,2, ... . (6.7) 
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The space $ is the collection of square integrable functions that belong to the maximal 
invariant subspace of H and for which the estimates ()6.7|) are finite. In particular, 
because ^^{r) satisfies the estimates ()(i.7j) . '^^{r) falls off at infinity faster than any 
polynomial of r: 

lim(l + r)>+(r) = 0, n = 0,l,2,.... (6.8) 

1 — >oo 

Obviously, the space $ can also be seen as the maximal invariant subspace of the algebra 
generated by the Hamiltonian and the operator multiplication by r. The estimates fl6.7j) 
are norms (see Proposition 1 at the end of this section), and therefore they define a 
topology (i.e., a meaning of convergence of sequences) r$ on 

ifl-^if-^ iff W'pI - V'^\\n,m ^0, n,m = 0, 1,2, . . . . (6.9) 

Once we have constructed the space we can construct its topological dual 
as the space of ^-continuous antilinear functionals on and therewith the RHS 
corresponding to the "in" states, 

$ C L2([0,oo),dr) C (6.10) 

One can show that \E^) indeed belongs to see Proposition 1 below. 

The conditions ()6.5|) and ()6.7|) that determine the space $ are very similar to the 
conditions satisfied by the Schwartz space on the positive real line, the major difference 
being that the derivatives of the elements of $ vanish at r = 0, a, h. This is why we 
shall write 

$ = 5(M+-{a,6}) , (6.11) 

where = [0, oo). With this notation, the RHS ()6.10p can be written as 

5(M+-{a, h}) C L2(R+, dr) C 5^ (M+-{a, h}) . (6.12) 

The "in" kets should be eigenvectors of H as in Eq. p.3|) . However, since the 
Hamiltonian acts in principle only on its Hilbert space domain, and since the "in" kets 
belong to the antidual rather than to the Hilbert space, we need to extend the action of 
H from $ to in order to specify how H acts on \E^). The theory of distributions 
provides us with a precise prescription for such extension: Given an operator A, the 
action of A on a functional |F) G is defined as 

{ip\A\F) = {A^^\F) , for all in $ . (6.13) 

It is important to realize that definition ()6.13|) makes sense only when $ is invariant 
under A\ 

A^^d^. (6.14) 

Note that when A is self-adjoint (e.g., A = H), then A^ = A, and when A is unitary 
(e.g., A = [/+), then A'f = A'K Definition (jOT^ can in turn be used to define the 
notion of eigenket of a self-adjoint observable: A functional |a) in is an eigenket of 
A with eigenvalue a if 

{(p\A\a) = {A(p\a) = a{{p\a) , for all in $ . (6.15) 
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When the "left sandwiching" of this equation with the elements of $ is understood and 
therefore omitted, we shall simply write 

A\a) = a\a). (6.16) 

Thus, within the RHS setting, the eigenvalue equation ()1.3p is to be understood as 

{ip+\H\E+) = E{^+\E+) , (6.17) 

This eigenvalue equation is proved in Proposition 1. 

By means of the unitary operator U^, which was constructed in Section I3?n we can 
obtain the energy representation of the space 

= (6.18) 

We shall denote the elements of by if~^{E) = ([/+(/?+)(£'), where E G [0, oo). Using 
the notation of Eq. ()6.1ip . the space will be also denoted as 

$+ = S+iR+-{a,b}), (6.19) 

and the energy representations of the triplets ()6.10|) and ()6.12|) will be respectively 
denoted by 

C L\[0, oo), dE) C (6.20) 

and by 

5+(M+-{a,6}) C L\R+,dE) C 5j(M+-{a,6}) . (6.21) 

The antidual extension of f/+ yields the energy representation of the "in" kets as 

\E+) = U+\E+). (6.22) 

It can be shown that {E'^) acts as the antilinear Schwartz delta functional, see 
Proposition 1. 

6.2. The "in" bras 

We now construct the bras that correspond to the kets \E~^). Likewise the definition 
of a ket, the definition of a bra is borrowed from the theory of distributions [21]. Given 
a function f{x) and a space of test functions the linear functional F generated by 
the function f{x) is an integral operator whose kernel is the complex conjugate of f{x): 



F{ip) = J dxip{x)f{x) . (6.23) 

By using prescription (j6.23p . we define the bra {^El associated with the "in" 
eigenfunction x+(r;i?) as 



{-^E\ip-^)= drx+(r;E)<^+(r), (6.24a) 
Jo 

which in Dirac's notation becomes 

POO 

{+E\^+)=l dr (+E|r)(r|v?+) . (6.246) 
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Here (+i?|r) denotes the complex conjugate of x'^i'f'] E). 

From definitions ()6.2aj) and ()6.24aj) . it follows that the action of the bras {'^E\ is 
complex conjugated to the action of the kets \E^): 



= {ip+\E+) = ^{E) . (6.25) 

In its turn, Eq. ()6.25|) shows that definition ()6.24aj) is meaningful when y?"*" G If we 
denote by the space of continuous linear functionals over $, then it is also easy to 
prove that {^E\ belongs to see Proposition 1 below. Therefore, the triplet 

# C L2([0,cx)),dr) C (6.26) 

is suitable to accommodate the "in" bras. Using the notation introduced in Eq. ()6.11|1 . 
we shall denote the triplet ()6.26p as 

5(M+-{a, b}) C L2(R+, dr) C 5'(M-{a, b}) . (6.27) 

Our next task is showing that the "in" bras are left eigenvectors of the Hamiltonian. 
For this purpose, we need to specify how H acts on the bras, that is, how H acts on the 
dual space The action to the left of an operator A on a linear functional {F\ G 
is defined as 

{F\A\^) = (F|AV) , for all (y9 in $ . (6.28) 

In turn, Eq. ()6.28|) can be used to define the notion of eigenbra of a self-adjoint operator: 
A functional (a| in $' is an eigenbra of A with eigenvalue a if 

{a\A\Lp) = {a\ALp) = a{a\Lp) , for all in $ . (6.29) 

When the "right sandwiching" of this equation with the elements of $ is understood 
and therefore omitted, we shall simply write 

{a\A = a{a\. (6.30) 

The eigenbra equation ()1.5|1 is therefore to be understood as 

{+E\H\ip+) = E{+E\ip+) , ip+ G 5(M+-{a, b}) . (6.31) 

Proposition 1 below proves this equation. 

The operator f/+ can also be extended to the dual space, and such extension can 
be used to obtain the energy representation of ('''-E'l, 

{+E\ = {+E\U+. (6.32) 

In Proposition 1, we show that {'^E\ acts as the linear Schwartz delta functional. 

It should be emphasized that there is a one-to-one correspondence between "in" 
bras and "in" kets, since to each "in" ket \E~^) there corresponds an "in" bra {'^E\, and 
vice versa. 
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6.3. The "out" kets 

The construction of the "out" kets closely parallels the construction of the "in" kets. 
By way of we define the "out" kets \E^) for each E G [0, oo) as 

\E-) : $ I — > C (6 33 ) 

^jj- I — > {^Ij-\E-) = dr^lj-{r)x'ir;E) , 

where x l'"! E) is the "out" eigenfunction of Eq. ()3.11|) . In Dirac's notation, the action 
of \E~) reads as 



{^Ij'\E-)= dr {ip-\r){r\E~) . (6.33&) 
Jo 

Similarly to the "in" case, the space $ in definition (jfi .330,11 must satisfy the 
following conditions: 

• The space $ is invariant under the action of H. (6.34a) 

• The elements of $ are such that the integral in Eq. (I6.33ap makes sense. (6.346) 

Likewise the x^(^) E), the E) behave like purely imaginary exponentials when E G 
[0, oo), which implies that conditions ()6.34a|) - fl6.346|) are equivalent to conditions ()6.3a|) - 
fl6.3&|l . Therefore, the space of "out" wave functions is the same as the space of "in" 
wave functions: 

G L2([0,oo),dr) I G P , < oo , n, m = 0, 1, 2, . . .} , (6.35) 

where the || \\n,m, are given by Eq. (j6.7p . It can be proved that the "out" kets \E~) 
belong to iS^(]R+-{a, 6}), and that they are eigenvectors of H, see Proposition 1 below. 

A comment on notation is in order here. We have used two different symbols, (f~^ 
and ip~, to denote the elements of one and the same space see Eqs. ()6.35|) and ()6.6|) . 
The reason why we use two different symbols is that we need to specify what kets are 
acting on the elements of When the elements of $ are acted upon by \E~^) {\E~)), we 
shall use the notation {(f^\E^) {(ip^lE^)). Another reason why we need two symbols is 
that, as we shall explain in Sec. [7| the time evolution of (f~^ is interpreted in a different 
way to the time evolution of ip' . 

By means of the unitary operator f/_, which was constructed in Section we can 
obtain another energy representation of the space 

f/_$ = $ . (6.36) 

We shall denote the elements of as ip~{E) = {U-il)~){E), where E G [0, oo). Using 
the notation of Eq. 1)6.111) . the space will be also denoted as 

= 5_(M+-{a,6}), (6.37) 

and the energy representations of the triples ()6.10|) and ()6.12|) . when obtained through 
[/_, will be respectively denoted by 

% C L2([0, oo), dE) c (6.38) 
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and by 

5l(M+-{a,6}) C L\R-^,dE) C 5!: (M+-{a, 6}) . (6.39) 
As well, the dual extension of f/_ yields the energy representation of the "out" kets, 

\E-) = U^\E-). (6.40) 
In Proposition 1, we prove that \E~) acts as the antilinear Schwartz delta functional. 

6.4- The "out" bras 

By using prescription ()6.23|) . we define the bra ("£"1 as 



{-E\i)~)= drx-{r]E)i)-{r) , (6.41a) 
Jo 

which in Dirac's notation becomes 

{-E\'ip^)= I dr {-E\r){r\il)-) , (6.416) 
Jo 

where {~E\r) denotes the complex conjugate of E). 

From definitions ()6.33a|) and ()6.41a|) . it follows that the action of the bras {~E\ is 
complex conjugated to the action of the kets \E~): 



{-E\i)-) = {i)-\E-) = i)-{E) . (6.42) 

In turn, Eq. ()6.42j) shows that definition ()6.41ajl is meaningful when G Therefore, 
the "in" bras belong to iS'(]R-{a, 6}), and they are left eigenvectors of H. Also, 
the energy representation of {^E\, = {^E\U-, acts as the linear Schwartz delta 

functional. 

The following proposition summarizes the results of this section: 



Proposition 1. The triplets of spaces \6.21 ) and llfi.l^} are rigged Hilbert spaces, and 



they satisfy all the requirements needed to accommodate the Lippmann-Schwinger bras 
and kets. More specifically, 

(i) The II ll^^m are norms. 

(a) The space S{M.~^-{a,b}) is dense m L^([0, oo), dr). 

(Hi) The space S(R~^-{a, b}) is invariant under the action of the Hamiltonian, and H is 
continuous. 

(iv) The kets |-E^) are antilinear functional over S(M.'^-{a,b}) , i.e., {E"^) G 
5^(M+-{a,6}). 

The bras {^E\ are linear functional over S(R~^-{a,b}) , i.e., (^-E| G 5'(M-{a,6}). 
(v) The kets {E"^) are (right) eigenvectors of H with eigenvalue E, 

H\E^) = E\E^). (6.43) 
The bras {^E\ are (left) eigenvectors of H with eigenvalue E, 

{^E\H = E{^E\ . (6.44) 
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(vi) In the energy representation, the "in" and "out" kets act as the antilinear Schwartz 
delta functional: 



(^|E+) = ^+(i?), (6.45) 



{^-\E-) =^-{E) , (6.46) 

whereas the "in" and "out" bras act as the linear Schwartz delta functional: 

= ^(E) , (6.47) 

{-E\i^~) =^'{E) . (6.48) 



The proof of this proposition can be found in Appendix C 



7. The time evolution of the Lippmann-Schwinger bras and kets 

In the previous sections, we obtained the (time- independent) solutions to the Lippmann- 
Schwinger equations. In this section, we obtain the time evolution of \E^), {^E\, ip~^ 
and ip~. 

7.1. The time evolution of the "in" states and kets 

In Quantum Mechanics, time evolution follows from the action of the operator e~'^*/^. 
This operator is unitary for each t G {—oo, oo), and the set 

|e-i^*/^| -oo<t<oo} (7.1) 

is a one-parameter unitary group. For each instant t, the time evolution of the "in" 
states is then given by 

<^+(r;t) =e-'^*/V(^), teR. (7.2) 

This time evolution can be conveniently written in terms of the "in" eigenfunctions by 
formally applying e~'^*/^ to both sides of the expansion (|5.18p : 



/■oo 

(e-^*/V) (r) = / dEe-'^*/ V(i^)x+(r; E) , 

^0 



(7.3) 



which in Dirac's notation becomes 



r\e-'"'/^\y^+)= I dEe-'^'*/'^(r|E+)(+E|<^+). (7.4) 



POO 

) = / dEe-'^'^^{r\E+){-^E\ip+) 
Jo 



Equation ()7.3j) is equivalent to defining e '^*/^ as the operator that in the energy 
representation acts as multiplication by 

(f/+e-'^*/V+) {E) = (e-'^*/V) {E) = e-^^*/V(^) • (7.5) 
Expressions ()7.3p and ()7.5p can be rigorously justified by way of Eq. ()A.5|) of Theorem 2. 
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In order to obtain the time evolution of the "in" kets \E^)^ we need to extend 
^-iHt/h antidual space. Such extension follows from prescription ()6.13|) : 

((/p+|e-'^*/^|E+) := (e'^*/ Vl^"") , '^^ e 5(R+-{a, h]) . (7.6) 

As noted in SeclHl this definition makes sense only when 4S(M+-{a, h}) is invariant under 

e'^*/^5(M+-{a, h}) C 5(M+-{a, h]) . (7.7) 

Such invariance is guaranteed for all t by Hunziker's theorem [22], see the theorem below, 
and therefore the time evolution of the "in" kets is well defined for all t. We note in 
passing that, as explained in Ref. [23], the invariance of the space of test functions is 
equivalent to having a well-defined Heisenberg picture, and therefore, from a physical 
point of view, it is clear that such invariance must hold. 

In order to state Hunziker's theorem, we need some definitions: A potential V is 
said to satisfy the Kato condition [24] if V{Hq) C V{V), and if there exist constants 
a < 1, 6 < cxD, such that, for all / G V{Ho), 

||\//|| <a||/7o||+&||/||. (7.8) 

When V satisfies ()7.8|) . then V can be seen as a small perturbation to the kinetic energy. 
For any positive integer n, we define a linear subset Dn of L^(M"'", dr) and a norm || ||„ 
on Dn by 

Dn:= fl Vir'^H"^), := sup \\r'H"^f\\, (7.9) 

, ^ k<n 

k<n ^— , 

^— I m<n—k 

m<n—k — 

where m is a positive integer. 

Theorem. (Hunziker) When the potential V satisfies the Kato condition ^7.^ , then the 
following holds for any positive integer n: 

(i) Dn is invariant under the unitary group e~'^*/^. 

(ii) For any f G Dn, Q~^^^l^f is continuous in t in the sense of the norm \\ ||„, and 
there exists a constant Cn such that 

||e-'^*/Vl|n<c„(l + |t|/n)"||/|U. (7.10) 

(Hi) For any f G Dn, 

the commutator being defined as 

\H, r"] = --r"P - — f" , (7.12) 

where r" and denote the first and second derivatives of r"^ with respect to r, and 
where P = —ihd/dr. In the L'^-norm, the integrand is continuous in r and bounded 
by constant (1 + |r|//l)""^' 
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We note that Hunziker's theorem, as stated in Ref. [22], is only vahd for 
Hamiltonians defined on L^(R^,dx) for any dimension A^. We therefore have to adapt 
the proof of Hunziker's theorem to our case, which involves a Hamiltonian defined on 
L^(]R"'", dr). Such adaptation will not be reproduced here, since it is straightforward. 

Our rectangular barrier potential clearly satisfies the Kato condition. Hence, 
Hunziker's theorem applies to our case. Because iS(M"'"-{a, 6}) = fl^o-^"' ^"^^ 
because Hunziker's theorem guarantees that each D„ is invariant under e~'^*/^, the 
invariance ()7.7|1 holds, and hence definition (j7.fi |) makes sense. 

After seeing that it makes sense, we are going to see that definition ()7.6|) yields the 
expected time evolution of the "in" kets: 

e-'^*/^|E+) = e-'^*/'^|E+) , -oo < t < oo , (7.13) 

which in the RHS language is to be understood as 

((^+|e-'^*/^|E+) = e-^^*/^((^+|E+) , (^+ G 5(M+-{a, h]) . (7.14) 

Equation ()7.14|1 follows from the following chain of equalities: 

((^+|e-^^*/^|E+) = {e'^'/^'<^+\E+) by definition ^ 

= e^^/^^+{E) by Eq. 

= e-^^*/^^+(^ by Eq. ^ 

= e-'^'/^{ip+\E+) . (7.15) 

Note that in the Hardy-function approach to the Lippmann-Schwinger equation, 
the time evolution of the Lippmann-Schwinger bras and kets is not defined for all times, 
but only for positive (or negative) times. Also, that the time evolution (j7.13j) is valid 
for all times could have been anticipated from the physics of a scattering process: The 
"in" solution of the Lippmann-Schwinger equation represents a monoenergetic ingoing 
particle prepared in the distant past (t = — oo) that hits the target and evolves into 
an outgoing particle in the distant future (t = +oo). Thus, physically, the process 
described by the "in" ket lasts from t = — oo until t = +oo, in agreement with ()7.13j) 
but in disagreement with the Hardy-function approach. 

7.2. The time evolution of the "in" bras 

The time evolution of {'^E] can be obtained by extending q-^^^/^ to the dual space. 
Such extension follows from definition ()6.28|) : 

(+^|e-i^*/^|(^+) := (+^|e^^*/ V) , e 5(R+-{a, h}) . (7.16) 

Likewise definition (j7.(i|l . definition ()7.1(ij) makes sense because 5(M+-{a, h}) is invariant 
under e'^*/^. 

By using Eqs. ()6.25|) and ()7.14|1 . it can be easily seen that definition ()7.16|) yields 

(+^|e-i^*/^|y,+) = e'^*/^+E|¥;+) , -oo < t < oo , (7.17) 
which, after omitting the v?"*", becomes the expected result: 

(+E|e-'^*/^ = e^^*/^(+E|, -oo<t<oo. (7.18) 
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7.3. The time evolution of the "out" states, bras and kets 

The time evolution of the "out" states can be obtained by formally applying e~'^*/^ to 
both sides of Eq. (fO^ : 

POO 

[e-'H'/^^-) {r) = dEe"'^*/^V^-(E)x"(r;E), (7.19) 
Jo 

which in Dirac's notation becomes 

noo 

{r\e-'^''^\ij-) = dEe-"^'/^{r\E-){-E\i,-) . (7.20) 
Jo 

Likewise Eq. ()7.Hj) . Eq. ()7.19j) is tantamount to defining e~'^*/^ as the operator that in 
the energy representation acts as multiplication by e"'^*/'^: 

([/.e-'^*/^-) {E) = (e-'^''^^^^ (E) = e-'^''^^-{E) . (7.21) 

Likewise expressions ()7.Hp and ()7.5j] . expressions ()7.19|1 and ()7.21|1 can be rigorously 
justified by way of Eq. ()A.5|) of Theorem 2. 

By using Eq. ()7.19|) . and following the same steps as for the "in" bras and kets, one 
can easily obtain the time evolution of the "out" bras and kets: 

(-E|e-'^*/^ = e'^*/^(-E| , -oo<t<oo, (7.22) 
^-iHt/h^^-^ = e-'^*/^|E-) , -cx)< t < oo . (7.23) 

As already mentioned in Sec. IHl the time evolution of the "in" states has a different 
interpretation from the time evolution of the "out" states. By the stationary-phase 
method, one can easily see that the "in" states are determined by the initial (i.e., 
prepared) condition that, as t — oo and r ^ oo, they move toward the potential 
region, whereas the "out" states are determined by the final (i.e., detected) condition 
that, as t ^ oo and r oo, they move away from the potential region (see also Ref. [25], 
p. 356). This is why we have denoted the "in" and the "out" states by two different 
symbols, even though they belong to one and the same space 

7.4. The time evolution of the free bras and kets 

In Ref. [17], we constructed the free bras and kets, but we did not obtain their time 
evolution. We do so here. 

Either by direct calculation, or by making the potential zero in the time evolution 
of the Lippmann-Schwinger bras and kets, one can easily see that the free bras and kets 
evolve in time in the expected way: 

{E\ e-'"°'/^ = e'^*/^(E| , -00 < t < 00 , (7.24) 
g-iHot/?i|^) = e-''Et/^\E) , -00 < t < 00 , (7.25) 
where the free bras and kets act as the following integral operators (see Ref. [17]): 



{E\^)= / dr^{r)xo{r;E), (7.26) 



{ip\E)= / drip{r)xo{r;E). (7.27) 
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8. Other results of scattering theory 

The "in" bras and kets can be used to expand the "in" states ip~^. This expansion is the 
restriction of the eigenfunction expansion ()5.18|) to the space # = 5(M"'"-{a, b}), 



{r\ip+)= dE{r\E+){+E\ip+) . (8.1) 
Similarly, by restricting Eq. ()5.15|) to $ = iS(M+-{a, 6}), we obtain 

/■oo 

= / dr(+E|r)(r|^+). (8.2) 
Jo 

The corresponding expansions of the by the "out" bras and kets follow from the 
restriction of the eigenfunction expansions ()5.33|) and ()5.32|) to iS(M+-{a, 6}): 

/■oo 

{r\^-)= dE{r\E-){-E\^-) , (8.3) 
Jo 

POO 

{-E\^-)= / dr {-E\r){r\^-) . (8.4) 
Jo 

Expansions dHU) and (Q are the way through which the RHS gives meaning to the 
formal expansions (j2.7|) and (j2.8|) . 

The S'-matrix element ()2.2j) can be written in terms of the action of the Lippmann- 
Schwinger bras and kets as in Eq. (j2.12j) . The expansion (|2.12|) plays an important role 
in resonance theory and is proved in [Appendix C By a similar argument to that used 
to prove Eq. ()2.12|) . one can also prove Eq. ()2.15|) . Many formal identities follow from 
Eqs. flTT^ and ^IJ^. For instance, 

{E\S\E') = {-E\E'+) = S{E)6{E - E') , (8.5) 

noo POO 

/ dr {E\r){r\E') = dr {^E\r){r\E'^) = 6{E - E') , (8.6) 
Jo Jo 

POD 

/ dr {~E\r) {r\E'+) = S{E)6{E - E') . (8.7) 
Jo 

As always, these expressions are to be understood within the RHS setting as part of a 
"sandwich" with well-behaved wave functions. 

The Lippmann-Schwinger equations and ()1.4j) are also understood as 

"sandwiched" with elements of For example, Eq. (jl.lll should be understood as 

= + (^I p ' (8-8) 

that is, 

POO POO 

dr^{r)x^{r]E)= I dr^{r)xo{r]E)+ l drds^{r)G'^{r,s]E)V{s)x^{s]E),{^.^) 



where if is an element of $ that can be attached a superscript + or — depending on 
whether it is an "in" or an "out" wave function. Note that since both \E^) and \E~) act 
on and since \E) is also well defined on the concerns raised in [10] do not appear 
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here. The concerns of [10] appear because in [10] it is assumed that the Lippmann- 
Schwinger kets are functionals over two distinct spaces of Hardy functions, whereas in 
this paper the "in" and "out" kets both act on one and the same space 

The M0ller operators Q± can be expressed in terms of the operators U± and Uq of 
Sec. El as [26] 

n± = uiUo. (8.10) 

As is well known, and as can be checked directly by using Eq. ()8.10|) . the M0ller operators 
intertwine the total and the free Hamiltonians: 

Ho = nlHn± . (8.11) 

Since our potential does not bind bound states, the M0ller operators are unitary 
operators on L^([0, oo), dr). The well-known expression for the ^-matrix operator in 
terms of the M0ller operators then reads as 

s = = u^U-ulUo . (8.12) 

The operator S is also a unitary operator on L'^{[0, oo), dr). In the energy representation, 
the operator ()8.12p acts as multiphcation by the function S{E) = J_{E)/J^{E). To 
be more precise, if we define the operator 5* as 

^:L2([0,oo),dE) ^ L2([0,oo),dE) 

/ ^ {Sf){E) = S{E)f{E), ^ ■ ' 

then it can be proved (see Appendix CD that 

S = UoSUo^. (8.14) 

The M0ller operators can be used to construct the space $o of asymptotic "in" 
and "out" V^""* states, 

$o = fi±*- (8.15) 
A vector 9?™ belongs to $0 if 

{+E\^^) = {E\^'-) , (8.16) 
where (p~^ = A vector belongs to $0 if 

{-E\^-) = , (8.17) 

where ip" = Q^ip"'^^. From the last two equations, it follows that 

n±\E) = \E^), (8.18) 

POO 

n± = / dE\E^){E\ . (8.19) 
Jo 

Again, Eqs. ()8.18|) and ()8.19|) are to be understood as part of a "sandwich" with elements 
of $ and $0- 
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9. Conclusions 

We have presented the RHS approach to the Lippmann-Schwinger equation. We have 
shown that the Lippmann-Schwinger bras are hnear functionals that belong to the 
dual space whereas the Lippmann-Schwinger kets are antilinear functionals that 
belong to the antidual space To every Lippmann-Schwinger ket there corresponds 
a Lippmann-Schwinger bra, and vice versa. The Lippmann-Schwinger bras (kets) are 
left (right) eigenvectors of the Hamiltonian, and their time evolution is defined for all 
times. 

The following diagram summarizes the results concerning the "in" states and kets: 
Ho; (p''^{r) $0 C L^([0, oo), dr) C #o 1-^) position repr. 

in+ in+ 

H; ip+{r) $ C L2([0,oo),dr) C \E+) position repr. (9.1) 

H; ^+{E) $+ C L'^([0,oo),dE) C |^+) energy repr. 

The results concerning the "out" states and kets are summarized by the following 
diagram: 



Ho; 


V'°"t(r) 




L2([0,oo),dr) 


C 


^0 


1^) 


position repr 


















H; 


ip~{r) 


* c 


L2([0,oo),dr) 


c 




\E-) 


position repr 


















H; 




c 


L\[0,oo),dE) 


c 




\E-) 


energy repr. 



Analogous diagrams summarize the results for the "in" and for the "out" bras. 

For the sake of simplicity, we have proved our results within the example of the 
spherical shell potential and for zero angular momentum. Nonetheless, with obvious 
modifications, the results of this paper remain valid for higher partial waves and for a 
large class of spherically symmetric potentials that includes, in particular, potentials 
of finite range. Therefore, we conclude that the natural mathematical setting for the 
solutions of the Lippmann-Schwinger equation is the rigged Hilbert space rather than 
just the Hilbert space. That rigged Hilbert space, however, seems to be unrelated to 
the rigged Hilbert spaces of Hardy functions of [8-14] . 
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Appendix A. The Sturm-Liouville theory 

In this paper, we have used several theorems that form the backbone of the Sturm- 
Liouville theory. These theorems can be found in the treatise of Dunford and 
Schwartz [20]. For the sake of completeness, we recall those theorems in this appendix. 

The following theorem provides the procedure to obtain the Green function of H 
(cf. Theorem XIII.3.16 of Ref. [20]): 



Theorem 1. Let H he the self-adjoint operator \4-4\l derived from the real formal 
differential operator \3.Ji\) by the imposition of the boundary condition \4-'Al - Let 
\m.{E) 7^ 0. Then there is exactly one solution x(r; E) of {h — E)a = square-integrable 
at and satisfying the boundary condition \4.i^ , and exactly one solution f{r;E) of 
{h — E)a = square-integrable at infinity. The resolvent {E — H)^^ is an integral 
operator whose kernel G{r, s; E) is given by 

r 2m x{r;E) f(s;E) ^ ^ ^ 

I IF WixJ) ' ^ f 5 
where W{x,f) is the Wronskian of x o-nd f 

W{x. f) = Xf - Xf ■ (A.2) 



The following theorem provides the operators U± (cf. Theorem XIII. 5. 13 of 
Ref. [20]): 



Theorem 2. (Weyl-Kodaira) Let h be the formally self-adjoint differential operator 
\3.Ji\l defined on the interval [0, oo). Let H be the self-adjoint operator Let A be an 

open interval of the real axis, and suppose that there is given a set {c'"i(r; E)., a2{r] E)} 
of functions, defined and continuous on (0, oo) x A, such that for each fixed E in A, 
{(Xi^r-jE), a2{r; E)} forms a basis for the space of solutions of ha = Ea. Then there 
exists a positive 2x2 matrix measure {pij} defined on A, such that 

(i) the limit 

{UfUE) ^ 



lim lim 



f{r)ai{r; E)dr 



(A.3) 
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exists in the topology of L'^{A, {pij}) for each f in L^([0, cxo), dr) and defines an 
isometric isomorphism U 0/ E(A)L^([0, oo), dr) onto L'^{A, {py}), where E(A) is 

the spectral projection associated with A; 

(a) for each Borel function G defined on the real line and vanishing outside A, 

UV{G{H)) = {[/,] e L\A, K}) I [Gf,\ e L'{A, {p,,})} (A.4) 

and 

{UG{H)f UE) = G{E){Uf UE), t = l,2,EeA,fe V{G{H)) . (A.5) 

The following theorem provides the inverses of U± (cf. Theorem XIII.5.14 of 
Ref. [20]): 

Theorem 3. (Weyl-Kodaira) Let H , A, {pij}, etc., he as in Theorem 2. Let Eq and 
El be the end points of A. Then 

(i) the inverse of the isometric isomorphism U o/£'(A)L^([0, 00), dr) onto L^(A, {pj^}) 
is given by the formula 

{U-'F){r)= hm hm H ( F^{E)a,{r■, E)p,,{dE)] (A.6) 

where F = [^1,^2] G L'^{A, {pij}), the limit existing in the topology of 
L2([0,oo),dr); 

(a) if G is a bounded Borel function vanishing outside a Borel set e whose closure is 
compact and contained in A, then G{H) has the representation 

poo 

G{H)f{r)^ / f{s)K{H,r,s)ds, (A.7) 

where 

2 r 

K{H,r,s) = Yl / G{E)ai{s;E)aj{r;E)pij{dE) . (A.8) 
i,j=i 

The spectral mcasmxs are provided by the following theorem (cf. Theorem XIII. 5. 18 
of Ref. [20]): 

Theorem 4. (Titchmarsh-Kodaira) Let A be an open interval of the real axis and O 
be an open set in the complex plane containing A. Let {ai{r\E)^ (T2(r;£')} he a set 
of functions which form a basis for the solutions of the equation ha = Ea, E E O, 
and which are continuous on (0, 00) x O and analytically dependent on E for E in O. 
Suppose that the kernel G{r, s; E) for the resolvent {E — H)~^ has a representation 



I E^- 1 d7-i{E)(Ti{r; E)aAs- E) , r <s, 

G{r,s;E)={ ^'''=' ''^ ' ^ ' (A.9) 
Y.h=,Otj{E)<J,{r-E)<j^{s-E), r>s, 
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for all E in Re(/7) (lO, and that {pij} is a positive matrix measure on A associated with 
H as in Theorem 2. Then the functions 9^ are analytic in Re(if ) fl O, and given any 
bounded open interval {Ei, E2) C A, we have for 1 < i, j < 2, 



liniA^n lim 



1 rE2~5 



I^:;sW,{E-ie)-9.^{E + ie)]dE 



lims^o lim,^o+ 2k lE^+si^tjiE - ie) - e±{E + ie)]dE . 



(A.IO) 



-e 
2 



k 

sin(fca) H cos(A;a^ 

IK 



Appendix B. List of auxiliary functions 

The coefficients in Eq. fj3.13|) are given by 
J,{E) = J,{k) 

J2{E) = J2{k) 

Js{E) = Js{k) 

ME) = J,{k) 
where k is given by Eq. ()3.10|) and k, is given by 



1 ■ / k 
-e""" I sin(A;a) - — cos{ka] 

2 \ iKj 



1 



1 . 



-ikb 



K 



K 



ikb 



1 - -^j e'^'Jiik) + (1 + ;^ ) e-'^'J2{k) 



K 



The coefficients in Eq. fl4.9jl are given by 



Al{E)=At{k) 



-e 
2 



k . 

1 ± - 1 e 



±ikb 



At{E)=At{k) = -e^^''[lT-]^ 



±ikb 



^f(E)=^f(fc) = -e-'^" 



1 

— ( 
2 

1 

-e 
2 



\ka 



At{E)=At{k) 
The coefficients in Eq. fl5.1?)|) are given by 

cos(A;a) 



1 + ^) e'^'^Atik) + (1 - ^) e-^'^'^^^ (A;) 
1 - ^) e'^'^Atik) + (1 + ^) e-'^'Mj (fc) 



cm 


= Ci{k) 


C2{E) 


= C2{k) 


C,{E) 


= Csik) 


cm 


= C,{k) 



sin(fca' 



2 \ ■ ' IK 

1 / k 

-e""" cos(A;a) + — sin(fca' 

2 \ IK 



-e 
2 



-ikb 



Jkb 



1 - -j e'^^'Cm + + ^ C2{k) 



(B.l 
(B.2 

(B.3 
(B.4 

(B.5 

(B.6 
(B.7 

(B.8 
(B.9 

(B.IO 
(B.ll 

(B.12 
(B.13 
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Appendix C. Proofs 

Here we prove some results we invoked tlirougliout tlie paper. In tlie proofs, it will be 
convenient to denote the dual and the antidual extensions of H (f/+) by, respectively, 
iU^) and H' {U'_^), in order to make clear that the operator H {U+) is acting outside 
the Hilbert space. 

Proof of Proposition 1. {i), (ii) and (m) were proved in Refs. [14,16]. 

(iv) From definition ()6.2aj) . it is easy to see that \E~^) is an antilinear functional. 
In order to show that {E'^) is ^-continuous, we define 

C+{E) := sup \x^{r;E)\ . (C.l) 

re[0,oo) 

From the expression for x^(r; E) in Eq. ()3.1ip . it is clear that C~^{E) is a finite number 
for each energy E. Now, because 



dr ip~^{r)x'^{r; E) 





oo 



< / dr|v?+(r)||x+(r;E)| 
Jo 

POO 

<C+{E) / dr|(^+(r)| 







C^{E) I dr-^(l + r)|v.+ (r)| 
1 + r 







r)2 



1/2 / /.oo \ 1/2 

dr 1(1+ r)Lp^{r) I 



= C+(B)I|43+IIlo, (C.2) 



the functional \E^) is continuous. In a similar way, one can prove that \E^) is also 
a ^-continuous antilinear functional. 

It is clear from definition ()6.24aj) that ("'"£'1 is a linear functional over iS(R"'"-{a, 6}). 
Because 

|(+^|^+)| = by (in2Sl) 

<C+(E)||<^+||i,o, bydlH (C.3) 

the functional {^E\ is ^-continuous. That {~ E\ is also a ^-continuous linear functional 
can be proved in a similar way. 

{v) In order to prove that \E^) is a (generalized) eigenvector of i/, we make use of 
the conditions satisfied by the elements of 5(]R+-{a, h}) at r = 0, oo: 

{ip+\H''\E+) = {Hip+\E+) 







(-^^ + ^W)^+(^)x^(^;^) 
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2m 



d(/?+(r) 



dr 



X+(r;E) 



+ 



oo 



+ / dry9+(r) 
'o 

r"00 





d2 



2m 



(^+(r) 



dx+(r;E) 



dr 



dr V5+(r) 



2m dr^ 

^2 d2 



+ V(r)U+(r;E) 



2m dr^ 

= E(v9+|i5+), (C.4) 

where we have used conditions and ()(i.8j) in the next to the last step and Eq. ()H.8j] 
in the last step. The proof that \E~) is an eigenvector of H follows the pattern of 
Eq. (in3I). 

On the other hand, because 



= {Hip+\E+) 



by (jFT^ 
by (Hn 

by (ra 



= E{ip+\E+) 

= E{+E\ip+) , by (C.5) 

the bra {'^E\ is a (generalized) left eigenvector of H. An argument similar to that in 
Eq. ()C.5j) shows that {^E\ is also a left eigenvector of H. 

{vi) Because 

= ^+(i?), (C.6) 

\E~^) acts as the antilinear Schwartz delta functional. A similar argument shows that 
\E~) also acts as the antilinear Schwartz delta functional. 
Equation ()6.47p follows from 

(+E|^) = (^+|E+) = ^(E). (C.7) 
Finally, Eq. ()(i.48j) follows from a chain of equalities similar to ()(I7j) . □ 



Proof of Eq. I{8.14\ )- order to prove Eq. ()8.14|1 . we first prove that 
{U^g){E) = S{E){U+g){E) , g e ^^([0, oo), dr) . 

Since by Eq. 

X^{r-E) = S{E)x~{r-E), 

and since 



(C.8) 
(C.9) 
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X~{r;E) = S{E)x^{r-E). (C.ll) 

By substituting Eq. (jdlljl into the integral expression ()5.32j) for the operator U-, we 
get to 



{U^9)iE) = / drg{r)S{E)x+ir;E) = SiE) dr g{r)xHr; E) . (C.12) 
Jo Jo 

Comparison of (|rn2)l with leads to (fTTs)) . 

Now, 

{UoSUo')f= {UoU^oU-UiUoUo')f by 
= {U-Ul)f 

= S{E){U+Ui)f hj^ 

= Sf, (C.13) 
which proves ()8.14j) . □ 

Proof of Eq. \2.1^) . Let ip'.ip'^ G Since ip' and v?"'" belong, in particular, to 

L^([0, oo), dr), we can let the unitary operator t/_ act on both of them. 

The vectors U-tp~ and U-^p^ belong to L^([0, oo), di?). Therefore, 

POO 

(V^-, ^+) = (f/_^-, f/_^+) = / dE (Z^r)CE)(f/_v5+)(i?) . (C.15) 

From Eq. (fTTsIl it follows that 

{U.v+){E) = S{E){U^V+){E). (C.16) 

Thus, 



(r,^+)=/ dE(^„^-)(E)5(i?)(f/+^+)(i?). (C.17) 
Jo 

Since ip' , (p~^ G we are allowed to write 



{U_ij-){E) = {^~\E-). (C.18) 
([/-,y.+)(E) = (+E|y.+). (C.19) 
Substitution of (fnT8|) and (fnT9|l into (fnTfj) leads to (jTT^ . □ 
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Figure 1. Schematic, pictorial representation of a scattering process. 



